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Abstract 

In this paper we define the concept of orthogonality between two factors 
"through another factor". Exploiting this property we have been able to ob- 
tain orthogonal main effect plans (OMEP) on non-orthogonal blocks requiring 
considerably smaller number of blocks than the existing methods. 

We have also constructed saturated partially orthogonal main effect plans 
(MEPs) for (i) an n 4 .2 3 experiment and (ii) an n 4 .2.3 experiment both on An 
runs. Here n is an integer > 3, n =/= 4. 

As particular cases, we have been able to accomodate four six-level factors 
on 8 blocks of size 4 each using the first method and on 24 runs using the second. 
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1 Introduction 



Main effect plans (MEPs) with mixed levels are often required for industrial 
experiments. The orthogonal main effect plans (OMEP) are, of course, the best 
option. However, due to the divisibility conditions, an OMEP for an asym- 
metrical experiments often require a large number of runs. For this reason, 
considerable attention has been devoted in the recent past to find main effect 
plans with the orthogonality condition relaxed to an extent. Such plans were 
first proposed in Wang and Wu (1992). Subsequently many others like Nguyen 
(1996), Ma, Fang and Liski (2000), Huang, Wu and Yen (2002) and Xu (2002) 
proposed and studied what they term "nearly orhogonal" arrays, concentrating 
mostly on two- or three-level factors. In many of these plans, there are factors 
non-orthogonal to three or more factors. As a result, in spite of the elegant 
combinatorial properties, the precision of the estimates go down. 

In our earlier paper [Bagchi (2006)] we have presented "inter-class orthogo- 
nal" MEPs, where a factor is non-orthogonal only to factors in its own class - to 
at most two factors in many of the plans. In this paper we continue the search 
for efficient MEPs with run size not too big. 

Study of an orthogonal main effect plan with possibly non-orthogonal block- 
ing was initiated in Mukerjee, Dey and Chatterjee (2001). They derived suf- 
ficient conditions for an OMEP to be universally optimal and also suggested 
a construction procedure for obtaining optimal OMEPs. In Bose and Bagchi 
(2007) we came up with a new set of suficient conditions for OMEP on non- 
orthogonal blocks of size two, requiring smaller number of blocks than Mukerjee, 
Dey and Chatterjee (2001). In this paper we generalise this idea. We define the 
concept of orthogonality between two factors "through another factor" (which 
may or may not be a blocking factor). Using this concept we have constructed 
an OMEP for a 4.n 3 experiment on n blocks of size 4 each, for n > 5. Treating 
the block factor as another treatment factor, one would obtain a 4.n 4 experi- 
ment on An runs. Although the main effects of this additional factor are not 
estimated with high precision, [see ( 3.14 )], there is considerable amount of re- 
duction in number of runs. [We may recall that an OMEP for an n 4 experiment 
requires at least n 2 runs]. The four-level factor may also be replaced by three 
two-level factors or one three and one two- level factor. 

We have presented another method showing that an OMEP for an s m ex- 
periment on bk non-orthogonal blocks of size k > 2 exists provided a connected 
binary block design with parameters (b, k, s) exist. Here m is the maximum 
number of constraints of an orthogonal array of strength two, k 2 runs and k 
symbols. In particular, we have obtained an OMEP for a 6 4 experiment on 8 
blocks of size 4 each. 
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2 An alternative to an existing MEP. 



We begin with a new main effect plan for a 3 4 2 3 experiment and compare it 
with the existing plan. Before that we need a few notation. 

Notation 2.1 A main effect (MEP) plan with a set F of m factors {A, B ■ ■ ■} 
and n runs will be represented by an array p(n, m; sa x sb • ■ ■), termed mixed 
array. In p rows will represent factors (in natural order) and Sa (respectively 
sa) will denote the set (respectively number) of levels of the factor A e F. 
Often, different rows have the same number of symbols and we may denote the 
plan as p(n, m; rii=i( s i) mi )> where 5^ i=1 m i — m - 

The vector of unknown effects of the levels of factor A will be denoted 
by the s^xl vector a. The replication vector of the factor A is the vector of 
replication numbers of its levels in the natural order and will be denoted by r A . 
R A will denote the diagonal matrix with entries as that ofr A in the same order. 
Ca will denote the coefficient matrix (C-matrix) in the reduced normal 
equation obtained by eliminating the effects of all factors other than A. 

For two factors A,B, the incidence matrix N A ' B is an sa x sb matrix 
with the (k,l)th entry as the number of occurences of the level k of A and I of 
B together. N A ' B = R A , when B = A. 

For a plan represented by a mixed array p(n, m; sa x sb • • ■) let Y denote 
the vector of yields. Then, our model is 



where p is the general effect, a is as stated above and Xa is a 0-1 matrix as 
described below. For A G F, the (i,j)the entry of Xa is 1 if the ith column of 
p contains j £ Sa in the row coresponding to factor A and otherwise. 

We present the well-known definition of an orthogonal array for the sake of 
completeness. 

Definition 2.1 An orthogonal array OA(n, m, si x • • • x s m , t), having m(> 2) 
rows, n columns, s\, . . . , s m (> 2) symbols and strength t(< m), is an m x n 
array, with elements in the ith row from a set of Si distinct symbols (1 < i < m), 
in which all t-tuples of symbols appear equally often (say A times) as columns 
in every t x n subarray. 

If Si = s, 1 < i < m, the OA is denoted by OA(n, m, s, t). A is called the index 
of the OA. Let us recall the nearly orthogonal array L' 12 (3 4 2 3 ) of Wang and Wu 

(1992), referred to here as AwwO-2)- 

We shall suggest an MEP Ai(12) for a 3 4 2 3 experiment described as follows. 




(2.1) 
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Let Ui = 



(2.2) 



Let V denote the array obtained by adding a row of all zeros to U\. Then, 



Ai(12) = 



V V + l V + 2 



[Here addition is modulo 3] . 

Comparing the MEPs AwwO-2) and Ai{12), we find the following. 

1. In both the plans, each two- level factor is orthogonal to every other factor. 

2. In both the plans, no pair of three-level factors satisfy proportional fre- 
quency condition. 

3. The C-matrices of the three-level factors for the plan Aww(12) are as 
follows. 



C Q = {7/3)K 3 , 



A,B,C,D. 



(2.3) 



[Here K n = I n — (1/n) J„, J is the matrix of all-ones.] 
Those for An {12) are as follows. 

C Q = 3K 3 , Q = A,B, C; Q D = {A/2,)K 3 . 

We observe that Ai(12) provides more information to factors A,B,C than 
Aww{^) but less to D. Further, C<3 1S bigger for A\{12), so that 

Q=A,B,C,D 

the total information is more for this plan. Thus, Ai(12) is useful in situations 
where all the factors are not equally important. 

Now, we go to the deeper question. We note that in A\yw{12), for each pair 
P i= Q,P,Q E {A, B, C, D}, Npq is the incidence matrix of a balanced block 
design (BBD), the "best" possible incidence matrix. On the other hand, for 
Ai(12), neither of Npq, P ^ Q is a "good" incidence matrix. How is then, Ca, 
for instance, is bigger for A\{12) ? 

Before going to the investigation of this mystry, we present two plans which 
are obtained by modifying A\{12) slightly. 



Let U 2 = 



1 
12 



and U 3 = [ 1 2 3 



(2.4) 



Then, 



M12) 



V 

u r 



V + l 

u,. 



V + 2 

Ui 



,i = 2,3. 



A 2 {12) and A 3 {12) are an MEPs for a 3 5 .2 and a 3 4 .4 experiment. The 
C-matrices of the first four (three-level) factors are, of course, same as those for 
Ai(12). The C-matrices of the new three-level factor and the two-level factor 
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of ^2(12) are 3^3 and 3K2 respectively. The C-matrix of the four-level factor 
in A3 (12) is 3K4. [K n is as in the statement below ( 2.3 ). 

To go into the mystry of "bigger" C-matrices of A,B,C in Ai{\2), we need 
some more notation. 

Notation 2.2 For L C F;U,V e F\L, C V y ;L will denote the (U, V)th subma- 
trix of the coefficient matrix ( C-matrix) in the reduced normal equation obtained 
by eliminating the effects of the member(s) of L. Cjj,U;L will be denoted by sim- 
ply Cjj-l- When L = S\U, Cu-l will be denoted simply by Cjj , to be consistent 
with Notation 2.1. 

Definition 2.2 If in a mixed array the factors A,B, and C satisfies the following 
condition, then we say that factors A and B are mutually orthogonal through 
C. 

N A,B = N A,C( R Cyl N C,B ( 3 5 j 

The following lemma is immediate from the definition above. 

Lemma 2.1 If in a mixed array a pair of factors A and B are orthogonal 
through C, then Ca,b,c = 0. 

Lemma 2.2 Suppose in a mixed array p two factors A and B satisfy the fol- 
lowing condition. 

For every Q E F = F \ {A, B}, A and Q are mutually orthogonal through 
B. Then, 

C A = C A , B . ( 2.6 ) 

Proof : We know that 

C A = C A -b - E a -b( h A;b)~(Ea;b) T , ( 2.7 ) 

where E A -,b = {{C A ,Q;B)) Qe p and h a-,b = {{Cp,Q;b)) p ^ q , PiQe p- 
From the hypothesis, Ea-.b = 0. Hence the result. □. 

Let us now look at the plans Aj(12),i = 1,2,3. We observe the following 
property. 

Lemma 2.3 For each of the plans ^(12), i = 1,2,3, the C-matrices of factors 
A,B,C satisfy 

Cp = C P ; D ,Pe{A,B,C}. (2.8) 
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Proof : It is easy to verify that each pair of factors (P, Q),P,Q 6 {A, B, C} 
, are mutually orthogonal through D. Hence the result follows from Lemma 2.2. 
□ 

Lemma 2.3 is the clue of the mystry of the "bigger" C-matrix in spite of 
the "poor" incidence matrices. We note that for A\yw(12) Ep.p ^ for each 
P 6 {A, B, C} and so, Cp < Cp.p>. As a result, even though Cp.p, is bigger for 
A WW (\T) than that for Ai(\2), Cp becomes smaller. 

However, no such facility is available for the factor D. Further, the incidence 
matrices of D with the other factors are not "good". This is why Cp, is smaller 
(than that in Aww{^)- 

Before going to the general construction in the next section, we present 
another example of a plan with two factors " orthogonal through another" . The 
plan is for a 3 3 experiment on 8 runs. 

~ 1 1 2 2 " 



A R = 



2 2 1 1 
20021001 
It is easy to verify that for As the condition obtained by interchanging A 
and C in ( 2.5 ) is satisfied, so that factors B and C are orthogonal through A. 
As a result, 



Cr — Cr-.a — 



-1 -1 
1 
3 1 



(2.9) 



We note that Cp has spectrum .1 .3 , which is only marginally "smaller" 
than the spectrum 1 .(9/4) 1 .3 1 of the C-matrix of a three-level factor with the 
hypothetical OMEP with best possible replication vector. 

Factor C also has the same C-matrix. 

Since A is non-orthogonal to two others, Ca is "smaller", as shown below. 



C A = (1/6) 
The spectrum of C A is O 1 .!^ 1 . 



4 -2 -2 
-2 7 -5 
-2 -5 7 



(2.10) 



3 Two series of orthogonal main effect plans on 
blocks of small size 

Notation 3.1 An MEP for an n"=i s i experiment laid out on b blocks of size 
k\.k,2, ■ ■ ■ kf, each will be represented by a mixed array p(n, m + 1; sa x sp, ■ ■ ■ x b; 
the last row of p represting the block factor. The set of factors will be represented 
by F U {bl}, bl denoting the block factor. R bl will denote the diagonal matrix 
with entries as k\.k2, ■ ■ ■ kb- 
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For the sake of easy reference, we write down the definition for orthogonality 
through the block factor for an equi-block-sized plan in a combinatorial form. 

Let us consider an MEP on b blocks of size k each blocks represented by 
a mixed array p(n,m + 1; sa x Sb • ■ ■ x b). Fix two factors A and B. For 
i 6 Sa,3 6 Sb 1 let Uij denote the number of columns in p with A at level i and 
B at level j. Let Bij denote the number of times the level combination (i,j) of 
A,B appear in the same block of p (need not be in the same column). 

Definition 3.1 Consider a mixed array p including a blocking factor of equal 
frequency k. Two factors A and B of p are said to be mutually orthogonal 
through the block factor if 

Bij = kUij, Vi£S A ,jeS B . (3.11) 

Note that Uij is nothing but the (i,j)th entry of JV A ' B , while Bij — 
N A,u N bi,B_ Thege together with the fact that R bi = kIb inl pi ies tnat ( 311 

) is a special case of ( 2.5 ) and hence Definition 3.1 is a special case of Defini- 
tion 2.2. 

Remark : The case k = 2 is considered in Theorem 3.1 of Bose and Bagchi 
(2007). 

We now generalise the arrays Ai(12),i = 1,2,3 to get an infinite series of 
MEPs. 

Theorem 3.1 (a) For an integer n > 5, there exist saturated MEPs for the 
following asymmetrical experiments on n blocks of size 4 each, 
(i) n 3 .2 3 , (ii) n 3 .2.3 and (Hi) n 3 A. 

(b) Plans (i) and (Hi) are orthogonal while (ii) is almost orthogonal. Specif- 
ically, in all these plans, the two, three and four-level factors are orthogonal to 
each of the n-level factors as well as the block factor. In plan (i) the two-level 
factors are orthogonal among each other. In plan (ii) the two and three-level 
factors are non-orthogonal among each other. 

(c) The two-level factors have C-matrix (fln)Ki in Plan (i) andnKi in Plan 

(ii) . The C-matrices of the three-level factor in Plan (ii) and the four-level factor 
in Plan (Hi) are given by nK^ and nK^ respectively. 

(d) The C-matrices of three n-level factors A,B and C are as given below. 

C Q = ((2 -1 ••• -1 )) , Q e G = {A,B,C}. 

Proof : (a) : Recall the arrays V and Ui, i = 1, 2, 3 described in ( 2.2 ) and 
( 2.4 . The following arrays represent the plans for the experiments (i), (ii) and 

(iii) respectively. Here the last row represent the block factor. [Recall Notation 
3.1]. 

M4n) = 



Ui Ui 
V V + l 



Ui 

V + n-l 



,2 = 1,2,3 
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(b) : It can be seen that each of the n-level factors satisfies proportional 
frequency condition with the two, three and four-level factors. 

(c) : These can be checked by straight computation. 

(d) : We note the following. Let the n-level factors be named as A,B,C in 
that order. Let G = {A, B, C}. Then, the incidence matrices between them are 
circulant matrices as given below. 

N PQ = L, P ^ Q, P,Q e G, (3.12) 

L is as given below. 



L = (( 2 1 ••• 1)) . ( 3.13 ) 

Further, one can verify that the pair of factors P and Q, P ^ Q, P,Q E G, 
are mutually orthogonal through the block factor. Now, it follows from Lemma 
2.1 that 



Cq = Cq-u,Q E G. 

Now, using ( 3.12 ) it is easy to verify that Cq, Q G G is as given in the 
statement. □ 

Remark 1: In the plans constructed in Theorem 3.1, the blocking factor 
may also be considered as a treatment factor named D say. Since D is non- 
orthogonal to each of A,B,C, it's efficiency would be low. Still it may be useful : 
particularly for the case n = 6 as using 24 runs we can accomodate four six-level 
factors, while on 36 runs one can accomodate at most three mutually orthogonal 
six-level factors. 

The C- matrix of D can be obtained as follows. 



C D =4I 4 -E D (H D )-(E D ) T , (3.14) 
where E D = [ M M M ] and Hd 



Ah L L 
L 4/ 4 L 
L L 4/4 

Here M is the circulat matrices described below and L is as in ( 3.13 ). 



M = 2(( 1 1 ••• 0)) ( 3.15 ) 

We now present another construction. 

Theorem 3.2 Suppose there exist a connected binary block design d with b 
blocks of size k each and v treatments. Let m denote the maximum number of 
constraints of an orthogonal array of strength two, k 2 runs, k symbols and index 
1. 

(a) Then, 3 an OMEP for a v" 1 ^ 1 experiment on bk blocks of size k each. 

(b) Further, the C-matrix Cp of every factor P is kCd- 
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Proof: Let us take the j th block of d consisting of the treatment set Tj 
say. Let Aj denote the OA(k 2 ,m, k, 2) with Tj as the set of symbols for the first 
m — 1 rows and the set {jk + 1, jk + 2, • • • , (j + l)k} for the last row. Then the 
array 



[ A, A 2 



A b 



is the mixed array p(bk 2 , m; v m ~ 1 .bk) representing the required MEP. 

It is easy to verify that condition ( 3.11 ) is satisfied by every pair of factor. 
Thus A represent an OMEP. Part (b) follows easily. 

Examples : 1. For a given v, there are many connected binary block designs 
with b blocks of size k each. Clearly, it is desirable that k is a prime or a prime 
power, so that one can accomodate k factors. Again, for the sake of economy, b 
should not be too large. Given these conditions, one possibility is to take 6 = 2 
and k the smallest prime power > [v/2). A binary connected design satisfying 
the above conditions always exist. Thus, Theorem 3.2 implies the existence of 
the OMEPs for the following experiments. 

(i) 4 3 and 5 3 on 6 blocks of size 3 each, (ii) 6 4 and 7 4 on 8 blocks of size 4 
each, (iii)8 5 and 9 5 on 10 blocks of size 5 each and so on. 

2. In the examples above the block design d is not equireplicate and hence 
the factors would not have equal frequency. To achieve equal frequency, of 
course, number of runs has to be bigger, particularly when v is prime. [ In that 
case the method of Mukerjee, Dey and Chatterjee (2001) might be useful]. We 
list the following connected and equireplicate block designs with composite v to 
be used in Theorem 3.2 together with the OMEP obtained. 

Block 1 Block 2 Block 3 



(a) v = 6 : b = 3, k = 4. d : 



1234 1256 3456 



OMEP : 6 4 experiment on 12 blocks of size 4 each. 



(b) v = 8 : b = k = 4. d: 



Block 1 Block 2 Block 3 Block 4 



1234 5678 1256 347 



OMEP : 8 4 experiment on 16 blocks of size 4 each. 

Block 1 Block 2 Block 3 Block 4 



(c) v = 10 : b = 4,fc = 5. d : 



12345 6789 10 12367 4589 10 



OMEP : 10 5 experiment on 20 blocks of size 5 each. 
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(d) v = 12 : b = 3,k = 8. d 



Block 1 



Block 2 



Block 3 



12345678 12349 10 11 12 56789 10 11 12 



OMEP : 12 8 experiment on 24 blocks of size 8 each. 

Remark 3 : It is interesting to find that the performance of the OMEPs 
derived from the block designs (a), (b), (c) and (d) above are quite good in the 
following sense. For each factor in each plan, the BLUEs of all the main effect 
contrsts except two have the same variance with the hypothetical (real in the 
case v = 8) OMEP without the blocking factor on the same number of runs. 
The varinces of the BLUEs of the two remaining contrasts are, of course bigger. 
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